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Abstract 

By assuming factorization of the GBD under the deconvolution integral for the hand- 



o 



bag diagram, we develop a method of solving this integral beyond the cross-over line. 
. As examples we use explicit models of deeply virtual Compton scattering (DVCS) 

£H | amplitudes to get solution for relevant GPDs. 
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1 Introduction 

Deeply virtual Compton scattering (DVCS) combines the features of inelastic processes 
^ ' with those of an elastic one. It has been realized [TJ[2j[3] that a straightforward generalization 
of the ordinary parton densities arises in exclusive two-photon processes in the Bjorken 
region, e.g. in Compton scattering with a highly virtual incoming photon, and in hard 
l/-) ! photoproduction of mesons. 

Generalized parton distributions (GPD) combine our knowledge about the one-dimen- 
sional parton distribution in the longitudinal momentum with the impact-parameter, or 
transverse distribution of matter in a hadron or nucleus [TJ |2j [3j H] . GPD cannot be measured 
directly, instead they appear as convolution integrals of the form 

m _ 

W, t ,V = [W^, (!) 
J x - £ + is 



where £ ~ — — — — is called skewness and x is the average longitudinal momentum fraction 

2 - x Bj 

of the struck parton in the initial and final states. This integral corresponds to deeply virtual 
Compton scattering (DVCS) in the "hand-bag" approximation, see Fig. [Ha), and x is the 
integration variable, not to be confused with the Bjorken variable XBj- 

Eq. (JTJ is an integral equation for the unknown function GPD. The solution is well 
known along the cross-over line x — £, where (see e.g. [I]) 

GPD(x = £, £, t, Q 2 ) = --Irm4(£, t, Q 2 ). (2) 

7T 

GPD (the nominator in Eq. is assumed to be universal for all exclusive processes 
of the type 7*p — > Vp, where V stands for a real photon or vector meson, the dependence 
on the processes coming from the "handle", calculable perturbatively. If so, many different 
processes of exclusive vector particle production can be described, according to Eq. ([TJ, by 
a universal GPD with different integration kernels, see Fig. [TJb). Since the GPD is not 
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Figure 1: a) Hand-bag diagram for DVCS; b) Hand-bag diagram for vector meson production 

known apriori, one seeks for models of GPD (see e.g. [5]) based on general constraints on 
its analytic and asymptotic behavior. The calculated scattering amplitudes (cross sections) 
are than compared with the data to confirm, modify or reject the chosen form of the GPD. 

In this paper we use two explicit models of DVCS and vector meson production (VMP) to 
obtain exact analytic expressions for GPD. To this end we propose a method of deconvolution 
of integral equations with a singular kernel of the Cauchy type to enable to go beyond the 
crossover trajectory and explore the hitherto unknown properties of the GPD. 

Regge pole models provide an adequate framework to describe high-energy, low t scat- 
tering phenomena. Being part of the 5* matrix theory, however, strictly speaking, they are 
valid only for the scattering of on-mass-shall particles. Still, the successful application of 
the Regge pole models in describing the off mass-shall HERA data opened the way to their 
use deep-inelastic scattering, DVCS and vector meson production at HERA. The appear- 
ance of the "forbidden" Q 2 dependence was circumvent by calling the Q 2 -dependent Regge 
trajectories "effective" ones. A particularly simple and efficient Regge pole model [7j with 
Q 2 -dependent residues (vertices) will be used below to guide our deconvolution procedure. 

In a more advanced Regge-pole model |8j, to be used below, the Regge trajectories and 
the residues do not depend on virtuality. Instead, the amplitude contains two (or more) 
Regge-pole terms, whose relative weight depends on Q 2 , mimicking the multi-pole nature 
of the so-called QCD Pomeron. Of great interest is the use of alternative models [6]. It is 
important to note that the Regge pole amplitudes are complex functions, their phase being 
fixed by the residue function. Since these model satisfy the basic properties of the theory, 
yet they fit the data, we shall use them to demonstrate the merits of our deconvolution 
procedure. 

The paper is organized as follows: in Section 2 (Simple DVCS amplitude) we introduce 
an explicit DVCS amplitude whose merit it simplicity, necessary in testing the capacity of 
our deconvolution procedure; in Section 3 (Deconvolution) we describe the mathematical 
basis for our method; in Section 4 (Factorization) we introduce the factorization hypothesis, 
essential in deriving GPD. In Sections 5 (Reggeometric model of DVCS) we apply our method 
to a more complicated model for the amplitude. 

2 A simple DVCS amplitude 

For the sake of clarity we start with a very simply but efficient model of the DVCS 
scattering amplitude [7j. This model accumulate in a compact way the main properties of 
the expected DVCS amplitude. Namely, it is Regge-behaved, has the required behavior in s, 
t and Q 2 , yet it fits the DVCS data measured at HERA by the Hi and ZEUS Collaborations. 
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The kinematics of those data are such that they correspond to diffract ive scattering, and 
consequently the t channel exchange is dominated by a single Pomeron trajectory. 
The explicit forms of the relevant scattering amplitude is [7J: 

A(s, t, Q 2 ) rp ^ p = -A e bia ^e b2l3{z \-ts/so) ait] = -A ^ bl+L ^ +b ^ z \ (3) 
where L = \n(—is/s ), 

a(t) = a(0) - «i ln(l - a 2 t), (4) 

and 

(5{z) = a(0) — «i ln(l — a 2 z), (5) 

where z = t — Q 2 is a new variable introduced in Ref. [TJ. 

From Eq. ([3]) we get the real and imaginary parts of the DVCS amplitude: 



lmA DVCS (x Bj ,t,Q 2 Q ) = sm(^)G(t,Ql)(^L-y (t \ (6) 
ReA DVCS (x Bj ,t,Ql) = - cos( ^-)G(t,Ql)(-^-Y {t \ (7) 



where 

G(t,Q 2 ) = e b(<*(t)+W,Q 2 )), ( 8 ) 

Skewness is defined in terms of the Bjorken variable x B j as £ ~ 2-% B or v,v, > XB i ~ T+l~ 

Since the Q 2 dependence in the model of Ref. [7J may not follow QCD evolution (ambigu- 
ous in DVCS), we keep it fixed at some value Q 2 that may be associated with the "frozen" 
QCD coupling constant. It has no effect on the deconvolution procedure in question. 

In the deconvolution procedure, Eq. ([I]), the variables t and Q 2 appear merely as "pa- 
rameters", therefore, for simplicity, we rewrite A DVC s as: 

A DV cs^t,Q 2 ) = -e-^B(t,Q 2 ) ( 1+1)°® , ( 9 ) 



where 

5(t,Q 2 ) = G(t,Q^)(g) Q(t) (10) 

is a real- valued function. 

With these ingredients we can now proceed to the deconvolution procedure. 



3 Deconvolution 

By substituting Eq. fl9]) in Eq. (CQ) we get (the dependence on the "parameters" is 
suppressed) : 

-i 
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where a < 1 and B is defined by Eq. (jlOl) . This is a singular integral equation with a 
Cauchy-type unknown kernel. By using the Sokhotski rule = p.v. (^) — Z7r5(z)), we 

obtain: 

1 

- e-^B (i±l ) ° = p.v. / ^^dx - ?7 rGPD(£, 0. (12) 

-1 

Since GPD is a real-valued function, we obtain its important property, called the cross- 
over trajectory condition: 

GPD(x = C, £, t, Q 2 ) = --ImA DVCS (C, t, Q 2 )- (13) 

7T 

At this moment we need to make approximation to proceed. Assume that GPD function 
could be factorized 

GPD(x,£) = /(£Mx), (14) 

we get 

""f«r(~J =p.v.y— d.-^K). (15) 

-1 

This is an integral equation for 4>(x) to be resolved. As we know from the Sokhotski-Plemelj 
theorem [9], there exist functions $ + (£) and $~(£) in the complex plain of the £ variable 
such that: $ + (£) is defined for Im(£) > and is analytic for Im(£) > 0; $~(£) is defined for 
Im(£) < and is analytic for Im(£) < 0, $ + (£), while are continuous on K (except fir 

a few points of the integrable discontinuity), and 

1 

$ + (0 + $-(0 = -P.v. I ^-Jx, £ G (-1, 1), (16) 



iri J x — £ 



-1 



<5> + (x) - $-(x) = <p(x), ie(-l,l), (17) 

$+(e) = $-(e), eeM\[-i,i]. (is) 

Here Eq. ffl8|) is a self-consistency condition for open-loop contours, see [9], Sec. 14.3-11 (for 
a more profound treatment of the problem see [ID]). This condition means (Eq. fflTj) ) that 
<p(x) = for i6K\ [—1, 1]. An alternative interpretation of this condition we will discussed 
below. 

By substituting Eqs. ([TBI) and (|17|) into Eq. (fT5l ) we arrive at the Riemann problem: 



®-(0 = -f^( 1 -^) • ief-l.l). (19) 



The analytic continuation of $~(£) for Im(£) < if very simple: $~(£) = |^j7gy (~£ 
From Eq. ( 1181) we can also find $ + (£) as an analytic continuation of $~(£): 

* + (Z) = -f^(^) a > M0>o. (20) 



2^/(0 V £ 

Suppose /(£) has no zeros or poles within the interval (—1,1). In this case we can write 
the solution even without any knowledge of the exact form of /(£)• The function $ + (£) has 
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an order a zero at the point —1 and an order a pole at the point 0. The cut from — 1 to 
along R and Eq. ([181) almost define on (-1, 1) \ {0}: 



$ + (0 



lira 

2 B 



_ i it a 
2 



£e (-1,0), 



£G (0,1). 



(21) 



(22) 




Figure 2: 

complex £ 



The function 
plane 



on the 



27ri/(o v e 

For a better understanding of this result, look at this 
problem from a different perspective: Eq. (118j) may be 
interpret as if there existed only one complex function 
on the £ plane — $, while $ + and $~ are values of this 
function for Im£ > and Im£ < 0, respectively. $ is a 
holomorphic function with branching points at —1 and 
0, where it has an order a zero and an order a pole. 
To determine its values unambiguously on all complex 
plane, we cut the plane from —1 to (see Fig. (j2J)). By 
encircling the point —1, initiating in the segment (—1; 0) 
(contour A in the figure), the phase of function $ will 
change by a factor e~ 2ma . However, if starting point will 
be in the segment (0; 1) (contour B on the figure), the 
phase of $ will not change since the factors coming from 
the zero and the pole will cancel; each other. 
From Eq. ( Tl7|) we have: 



<p{x) = (1 



-2-Kia 



17TQ , 

e 2 B 1 1 + x 



' 2mf{x) 
(p(x) = 0, 

Finally, for the GPD we obtain the solution: 



x 



, x e (-1,0), 
x e (0,1). 



GPD(s,0 = (1 - e" 2 -) e " o lf/i ( — ) -re (-1.0). 



2nif(x) \ x 



GPD(x,O = 0, 



x e (0,1). 



(23) 
(24) 

(25) 
(26) 



We see that this solution is unphysical, because GPD must be a real- valued function for 
physical values of £ (0 < £ < 1, see [3]). In the next section we discuss how to satisfy this 
condition. 



4 Factorization 

In the previous section we have used the factorization assumption for the deconvolution 
of Eq. ([1}). At this point it is hard to say how good is this assumption. We shall estimate it 
by its consequences. 

As we have seen in the previous section, the solution for GPD critically depends on the 
analytic properties of the left-hand side of the Eq. ( 1151) for £ G (—1; 1). Our degree of 
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freedom is the function /(£) in Eq. (j!4p . Its zeros and poles may help us in getting physical 
values for GPD. 

Let us construct a solution that will be a real- valued function for x G (a, 6), where (a, 6) 
is an arbitrary segment along the real axis. It is easy to check, that the function 

/(£)= f^Y(£-o)*(&-£)-fy0, (27) 



solves this problem. Here is an arbitrary real- valued function without zeros or discon- 
tinuities of fractional order along [—1,1]. Let us check this. 
By substituting ( )27|) into ( !T9|) . we get: 



*-(0 = -!p77T(£-ar f (&-o f , ee(-u). (28) 

2-Kigig) 

By repeating the procedure from the previous section, we see that the expression for $ + 
(g(£) has only zeros or discontinuities of integer order, so it changes phase only by integer 
number of 2m) : 

$ + (0 = -|^y(e-a)- f (&-O f , (29) 



lira 

$ + (0 = "|-^y(e-a)- f (&-O f , ^R\[a,b}. (30) 
From Eq. (ITTj) we have: 

</?(x) = — sin f — ") — — (x — a)~?(b — x) 2, 16 (a,6), (31) 
V 2 / 7T<7(:e) 

y?(x) = 0, i6R\[a,i)]. (32) 

Accordingly, for GPD we have: 

rpn/ t . • Bp(Q ^-a)(6-x) V fl + tY , , n 

GPD(o;, ~ sm (~2~j ^gjxj \ (x — a)(b - fl J \ £ / ' ^ )' ( 33 ) 

GPD(x,O = 0, ieR\[o,l], (34) 

The most interesting case for us is that ofa = — 1, 6 = 1, then GPD is given by the 
expression: 



GPD(x,£,t,Q 2 ) = - sin 



na(t)\ B(t,Q*)g(£) m+£)(l-z)\^ fl + p*® 



7T g(x)\(l + x)(l-0 



ng(x) V(l+x)(l-0 



a(t) 



ImA DVCS (Z,t,Q 2 ). (35) 



This GPD satisfies the cross-over trajectory condition ({TBI . It is easy to see that any 
other function / except ( )27|) will violate the reality of GPD. So, this is a general solution. 



6 



5 Reggeometric model of DVCS 



Another simple and explicit albeit more advanced model for DVCS, as well as for vector 
meson production (VMP) was proposed recently [8|. 

Here the Pomeron is also unique for all reactions, but the scattering amplitude contains 
two terms, a "soft" (s) and "hard" (h) one, "weighted" by Q 2 -dependent pre-factors: 



A(s,t,Q 2 ,M 2 v ) = , \ . na e-< W«) (—) asit) e 2 {^ + ^4 ) '+ 



+ _ r^_ e -*W*)(_f_) e W . (36) 



1 + « 



In a sense, it mimics the multi(infinite)-component QCD Pomeron and, in principle, is 
applicable to any exclusive reactions, be it "soft" or "hard". The parameters of the model 
are fixed from "first principles" and/or from the fits to the data, see [8]. 

Similarly to Sec. 3, and by using the relation s = ® ' we rewrite this amplitude in 
terms of the variable characteristic of GPD: 

t, Q\ Ml) = e-^B s (t, Q 2 , M 2 V ) (i±l) ^ + e -^B h (t, Q\ M 2 V ^ 1 + ^ ^ 

(37) 

where 



B l {t,Q\Ml)= , \ (^®A%^)\ i = s,h. (38) 



It is easy to see that our simple factorization approach is not working here any more. 
Really, let assume factorization to be of the form 

GPD(x, e, t, Q 2 , Ml) = f(£, t, Q 2 , Ml)(p(x, t, Q 2 , M 2 V ), 

with an arbitrary function /. The relevant solution will appear in the form (the dependence 
on the "parameters" t,Q 2 ,Ml is again hidden): 



GPDfoO = -f. ( (1 - e -^K^/(0)) e -^ Bs ( 1±^V S + 
2mj\x) \ ' \ x J 



+ (1 - c -™(«*-*vtt))) e ~^B h (^Y) 



(39) 



where «/(£) is an additive phase shift depending on the analytic properties of /. Let a s ^ ah 
and their values be independent. Then we can arrange the GPD to be a real- valued function 
for a certain point x = xo,t = xq, both terms in ( 139]) complex, their being real. However, 
by changing the value of x, we break this condition, since atf cannot depend on x. 

We may proceed in the following way: the equation for GPD (pQ) is linear, so we can split 
our amplitude (|3T|) in two parts and search for the solution as the sum of solutions for both 
parts. 
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Let t, Q\ Ml) = A a (£, t, Q\ M 2 V ) + A h {£, t, Q\ Ml), where: 



A s (U,Q 2 ,Ml) = e-^B s (t,Q\Ml) (l+lY 3it \ (40) 



a h (t) 

A h {t, t, Q 2 , Ml) = e-^B h (t, Q\ Ml) [ ^ ) . (41) 



By using the known solution ( 1351) for (0) we write: 



GPD,(,,{, . . . ) = sin (!2i) ^ + ~ fi±«r • 

v 2 ; i s ,(i) V(i + i)(i-oy v { / 

GPD k («, £,...)=«.(— )-_( ITT - yIr - Iy j , (43) 

where g s (£,,t,Q 2 , Ml) and gh(£,,t,Q 2 , Ml) are two independent arbitrary real-valued func- 
tions that have only integer order zeros and/or singularities for £ e [—1,1]. 
The resulting GPD will be the sum of ( gg} and f l43|) . 



6 Conclusions 

We have suggested a method for deconvoluting general parton distributions beyond the 
cross-over line based on the assumed factorization properties of the GPD under the convo- 
lution integral. We test the method by using two simple, explicit Regge-pole based models 
of the DVCS amplitude. It is important that our scattering amplitudes, due to the presence 
of the Regge phases, are complex-valued functions. Otherwise, the phase can be recovered 
either by dispersion relation methods or from the interference between the DVCS amplitude 
with Bethe-Heitler processes, see |1 1 j . We hope that this method can be applied to other 
models of DVCS (or vector meson production (VMP)). Comparison with alternative models 
of DVCS (VMP) and the resulting GPDs is of great interest. 

Although the models [7] and [8], used in the present paper, contain some Q 2 variation 
(introdeced phenomenologically) , its value was fixed for simplicity at some Qq that can be 
interpreted as a fixed QCD running constant. 
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